/AD-A185  487 


UNCLASSIFIED 


COMPARING  DISPERSION  EFFECTS  AT  VARIOUS  LEVELS  OF 
FACTORS  IN  FACTORIAL  EXPERIMENTS^)  CALIFORNIA  UNIV 
RIVERSIDE  DEPT  OF  STATISTICS  S  GHOSH  EE  AL.  AUG  87 
TR-159  AFOSR-TR-87-1086  4AFOSR-87-0048  F/O  12/2 


AD- A 185  407 


REPORT  DOCUMENTATION  PAGE 

•7 clasTI FICATION  MMt  lib.  RESTRICtTveTmA r'kings 


I’ll.  ntfOHT  SECURITY  CLASSIFICATION  MVWW 

’  I _ Unclassified _ I  J  |  | \ 

'  Jj».  SECURITY  CLASSIFICATION  AUTHQflUr  pi 

...  .  J  NA _ flELECTE 

t&fbiijpa  j*  2b.  DECLASS. F, CAT. ON, OOVYNGHAO.S»|HEg£f  0  6  1987 
I  NA  _  Yflfl  __________ 

■I  K  4.  PERFORMING  OFCANI2A  I  IQI I  III  I'j^^MuMSERIS) 

I  Technical  Report  No.  159 


“TOGfitErCOET 


S3.  DISTRIBUTION, AVAILABILITY  of  report 

Approved  for  public  reLics:-, 
_ distribution  unlimited _ _ 

5.  MONITORING  ORGANIZATION  REPORT  NUMBERISI 


1 _  U _ AFQSK«Tll«  87-106 

I's*.  NAME  OF  PERFORMING  ORGANIZATION  >b.  OFFICE  SYMBOL  7*.  NAME  OF  MONITORING  ORGANIZATION 
\  TIniuorci  fw  nf  (7a  1  i  f  ornia  Ilf  applicable)  a  rnnn  /vrw 


|  University  of  California  "f‘ 

jj  Riverside 

[*6e.  AOORESS  Idly.  Stott  and  ZIP  Codtl 

•'?!  Department  of  Statistics 
'  I  University  of  California,  Riverside 
rv'Vv  t  Riverside,  CA  92521 


AFOSR/NM 

Tb.  AOORESS  icily.  Slott  and  ZIP  Codtl 

Bldg.  410 
Bolling  AFB 
DC  20332-6448 


SO.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

AFOSR-87-0048 


Ft.  AOORESS  tCII y.  Slalt  aid  ZIP  Codtl 

31dg .  410 
Bolling  AFB 

DC  20332-6448  _ 


u.  title  tinciudt  Steamy  clarification)  Comparing  Dispersion 
effects  at  various  levels  of  factors  in 


,i  iz.  personal  authorisi  factorial  experiment, 
j  Subir  Ghosh  and  Eric  E.  Lagergren _ 

'■  5  13a  TYPE  OF  REPORT  1 3b.  TIME  CO  V  E  R  E  O 

I  Irrty rTlit  from  12/86  to  08/87 

1'  "  IB.TuFPLEMENTAR  Y  NOTATION 

Submitted  to  Technometrics 


|  10.  SOURCE  OF  FUNDING  NOS. 

PROGRAM 

PROJI  CT 

TASK 

WORK  UNIT 

element  no. 

NO. 

NO. 

NO. 

U 1CQF 

£304- 

_ 

1«.  DATE  OF  REPORT  (Yr..  Mo..  Day) 

August  1987 


15.  page  COUNT 

17 


I  _ _ _ — — — - - - ~ 

[  )7  COSATI  COOES  11  SUBJECT  TERMS  .Conllnu*  on  reverie  if  nectuary  and  identify  by  Moc*  number) 

\  field _ group _ sub.  gR. _  Adjusted  residuals.  Design,  Dispersion  effects,  Error, 

|  ~~ _ _  Factorial  experiments.  Linear  models.  Quality  control. 

I  1 _ __ _ _ _ — - - 

J  19.  ABS  TRACT  tConlinut  on  reverie  if  nectuary  and  identify  by  blocb  number) 

-  This  paper  is  an  attempt  to  understand,  measure  and  compare  dispersion  effects  at  dif¬ 
ferent  levels  of  factors  in  factorial  experiments.  The  simplest  setting  is  considered 
in  order  to  develop  better  comprehension  and  insight.  The  properties  of  the  proposed 
descriptive  measures  are  examined.  A  method  of ^"adjusting  residuals  and  its  use  in 
comparing  dispersion  effects  are  discussed.  Illustrative  examples  are  also  given.  The 
problem  considered  in  this  paper  arises  in  quality  control  studies  and  the  methodologies 
are  applicable  to  industrial  experiments. 


?0~OiSt  Rl  0UTION/A  VA»LA8>LlTY  of  ABSTRACT 

'  (NO.  ASSl*  1  E  O/UNUMi  TE  O  0  SAMg  AS  RPT  CJ  OTIC  USERS  CD 


87  9  24  214 

21  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


33b.  TELEPHONE  NUMBER 
tinciudt  1  rra  Code) 

rmi'i  7F,7-sn?7 


23c.  OF  F  ICE  SYMBOL 

AFnsR/NM 


5 


i  i  .  <  .  *  .  i 


i  ' 

f. 


«/  '  AS'  v-  •■  '  ;  •-  - 

,y,  4te$*$’nsk:x  k\. : ;  ?,•  • 

..-•v  '  r  /  \  •)  ■  •  if 


'  //•  i 


V;  ;  •' 

<  4 ./  ■:  : 

■•>•••  ■:  : 

’  '  '•  v? ,/  - 

;  -• 

:  \  ■  /  C«  \ 

'  I  <  .  s  ?  * 

'  n>0:ir 7:  : 

>■  I  'i.  •  v 


;  *  •  /  ■» 

.  f  .  : 


i  v  ;  -  , 

\  '  r  ■ .  • 


n'v 


•J  * 


AFOSR.tr.  8  7-  1086 


1 


I 


Comparing  Dispersion  Effects  At  Various  Levels 
Of  Factors  In  Factorial  Experiments 


Subir  Ghosh  and  Eric  S.  Lagergren 


Technical  Report  No.  159 


:  ...  1  ■  TVf-.-  1  f  1 


• «  t  • "  v*  t  1  V*  ;/*,  •  *  »V  \ 


v  f  !■  S  '  .  X 


*  -  •  •!  •’  s  ••  •' 


\  *  *•••'*  1 


j  •  .  .  •• ; -I  •  ;  ; ,  >)  ■'  '  v.'v 

:  .  -  /  -  j  ■  *  .  »  v  *  •  »  *  •  -*  !  *  -  > 


University  of  California 


Riverside  %' 


Department  of 
Statistics 


Comparing  Dispersion  Effects  At  Various  Levels 
Of  Factors  In  Factorial  Experiments 

Subir  Ghosh  and  Eric  S.  Lagergren 

Technical  Report  No.  159 


Department  of  Statistics 
University  of  California 
Riverside,  CA  92521 


August,  1987 


Accesion  for 


By . 

Di*'-tr  ib«:t*ort  / 


A-.UilO!  i::ty  CofitS 

|  A*. fill  .  or 
womi 


Copy 

vINSpCrrro 


Comparing  Dispersion  Effects  At  Various  Levels  Of  Factors  In 

Factorial  Experiments 


Subir  Ghosh*  and  Eric  S.  Lagergren 
University  of  California 
Riverside,  CA  92521 


Summary 

This  paper  is  an  attempt  to  understand,  measure  and  compare  disper¬ 
sion  effects  at  different  levels  of  factors  in  factorial  experiments. 

The  simplest  setting  is  considered  in  order  to  develop  better  compre¬ 
hension  and  insight.  The  properties  of  the  proposed  descriptive  measures 
are  examined.  A  method  of  "adjusting"  residuals  and  its  use  in  comparing 
dispersion  effects  are  discussed.  Illustrative  examples  are  also  given. 
The  problem  considered  in  this  paper  arises  in  quality  control  studies 
and  the  methodologies  are  applicable  to  industrial  experiments. 

Key  Words:  Adjusted  residuals,  Design,  Dispersion  effects,  Error, 
Factorial  experiments.  Linear  models,  Quality  control. 
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1.  Introduction 


An  Important  problem  In  quality  control  studies  is  to  find  an 
optimum  combination  of  levels  of  control  factors  in  achieving  stability 
against  noise  factors  (see  Taguchi  and  Wu  1985).  Both  "location"  and 
"dispersion’'  effects  of  levels  of  factors  are  pertinent  to  measure  from 
the  data  in  resolving  this  problem.  This  article  considers  the  problem 
of  measuring  dispersion  effects  of  levels  of  factors  in  both  replicated 
and  unreplicated  factorial  experiments.  The  concept  of  dispersion 
effects  in  factorial  experiments  was  introduced  in  Taguchi  and  Wu  (1985) 
for  replicated  factorial  experiments  and  in  Box  and  Meyer  (1986)  for 
unreplicated  factorial  experiments.  Factorial  experiments  may  be 
complete  or  fractional  factorial  under  completely  randomized  designs. 
Although  for  clarity  we  consider  2m  factorial  experiments  in  this 
article,  the  ideas  presented  can  easily  be  generalized  to  any  symmetric 
or  asymmetric  factorial  experiments.  Kackar  (1985),  Phadke  et.  al. 
(1983)  and  Nair  (1986)  made  pioneering  contributions  to  this  area  of 
research.  Ghosh  (1986)  used  the  search  linear  models  (see  Srivastava 
1975)  to  explain  dispersion  effects  in  factor  screening  experiments. 

We  first  assume  that  for  the  fitted  model  to  the  data  there  is  no 
significant  lack  of  fit.  We  then  propose  three  measures  of  dispersion 
effects  at  levels  0  and  1  of  m  factors.  All  three  of  them  are  relevant 
in  replicated  factorial  experiments  and  two  of  them  are  applicable  to 
unreplicated  factorial  experiments.  We  observe  that  the  measures  of 
dispersion  effects,  based  on  residuals  obtained  by  the  least  squares  fit 


of  Che  model  to  the  data,  at  levels  0  and  1  of  a  factor  are  correlated  In 


moat  situations.  We  introduce  a  method  of  adjusting  residuals  and  then 
propose  measures  based  on  residuals  and  adjusted  residuals. 

2.  Dispersion  Effects 

We  consider  a  2m  factorial  experiment  under  a  completely  randomized 
design.  Let  T(nxa)  be  the  design.  The  rows  of  T  denote  treatments  and 
the  columns  denote  factors.  The  design  T  is  called  an  inner  array  for  m 
controlled  factors.  For  various  level  combinations  of  noise  factors 
(outer  array),  we  get  replicated  observations  for  every  treatment  in  T 
(see  Taguchi  and  Wu  1985).  In  the  experiment,  we  take  r  (>  1)  observa¬ 
tions  for  every  treatment.  The  case  r  ■  1  is  called  the  unreplicated  ex¬ 
periment  and  the  case  r  >  1  is  called  the  replicated  experiment.  Again, 
for  simplicity  equal  replication  is  considered  for  the  replicated  experi¬ 
ment  and  the  idea  is  easily  extendable  to  unequal  replications.  Let  y 
be  the  jth  observation  for  the  ith  treatment,  y^  be  the  mean  of  all 
observations  for  the  treatment  i,  i*l,...,n  and  j»l,...r,  and  (N  ■*  nr)  be 
the  total  number  of  observations.  The  standard  linear  model  for  the 
experiment  is 


E(jO 

-  X8, 

(1) 

V(£) 

-  a2I, 

(2) 

Rank 

X  -  p, 

(3) 

where  jjr(Nxl)  is  the  vector  of  observations,  j5(p*l)  is  the  vector  of 

factorial  effects  considered  in  the  experiment,  X(Nxp)  is  a  known  matrix 

2 

that  depends  on  the  design  T  and  o  is  an  unknown  constant.  We  denote 
H  ■  X(X'X)  *X'  and  R  ■  (I-H).  The  vectors  and  are  the  vector  of 
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least  squares  fitted  values  and  the  vector  of  residuals,  respectively. 

The  fitted  values  for  all  observations  corresponding  to  the  ith  treatment 

A 

are  identical  and  is  denoted  by  y^,  i-1 . .  Suppose  that  for  the 

fitted  model  to  the  data  there  is  no  significant  lack  of  fit.  The  sum  of 

nr  * 

squares  of  error  is  SSE  ■  I  E  (y.,-y, )  »  the  mean  square  of  error  is 

i-1  j-1  J 

MSE  -  (SSE/(N-p)),  the  sum  of  squares  of  pure  error  is  SSPE  * 
nr 

E  E  (y  -y  )  and  the  mean  square  of  pure  error  is  MSPE  - 
i-1  j-1 

(SSPE/n(r-l)).  Note  that  both  MSE  and  MSPE  are  measures  of  error  vari- 
2 

ance  o  .  We  now  take  MSE  and  MSPE  as  descriptive  measures  of  noise.  We 
then  express  MSPE  as  the  weighted  average  of  (MSPE)^  and  (MSPE)^,  where 
(MSPE)u  is  called  the  contribution  of  the  level  u  (u  -  0,1)  of  the 
factor  to  MSPE.  Formal  expressions  of  (MSPE)U,  u  -  0,1  are  given  in  the 
next  section.  We  do  the  same  for  MSE.  Different  levels  of  a  factor  may 
contribute  differently  to  MSE  and  MSPE.  In  general  the  contributions  of 
levels  of  a  factor  to  noise  (measured  by  MSPE  or  MSE)  are  called  the 
dispersion  effects  of  levels  of  the  factor.  The  main  theme  of  this  paper 
is  to  investigate  the  possible  ways  of  measuring  and  comparing  dispersion 
effects  of  levels  of  factors.  We  would  like  to  make  it  clear  that  we  are 
not  presenting  a  population  model  for  dispersion  effects  and  the  proposed 
measures  in  this  paper  are  all  descriptive. 

3.  Measuring  Dispersion  Effects 
3.1.  First  and  Second  Sets  of  Measures 

We  take  a  single  factor  out  of  m  factors  and  develop  the  methods  of 
measuring  dispersion  effects  at  the  level  0  and  1  of  the  chosen  factor. 

We  do  not  introduce  any  notation  for  the  chosen  factor.  This  is  to  keep 
our  presentation  neat  and  clean.  We  define  for  i-l,...,n. 


-4- 


1  if  the  level  of  the  factor  in  the  ith  treatment  is  1, 


0  if  the  level  of  the  factor  in  the  ith  treatment  is  0. 


Let  (NxN)  be  a  diagonal  matrix  with  n  sets  of  diagonal  elements  and 


the  elements  in  the  ith  (i-l,...,n)  set  are  equal  to  6..  Ue  define 


Dq  »  I-Dj .  It  can  be  seen  that  D^Dq  -  0  and  both  and  Dq  are 


idempotent  matrices.  We  have 


SSE  «  j£'  R  D  R  x  +  i.'  RD.Rj 


n  r 


E,  *  +  E.  • 


i-1  j-1 


i-1  j-1 


n  r 


SSPE  -  z  z  «  (y  -y  )2  +  E  I  (l-0(y..-y.)2. 
i-1  j-1  2  i-1  j-1  1  1 

The  first  set  of  measures  of  dispersion  effects  of  levels  of  the  factor 


n  r  ,  -  ' 
E  E 

s?(D  •  j-1-  1 


(  E  O  (r-1) 

i-1 


s*u> 


E.  ,E. 


(  e  (c-n 

i-1 


at  the  levels  1  and  0,  respectively. 


We  have 


S2(l)  +  \  — 


Z  (1-6  ) 


s20(i) 


<■ 

*; 


$ 

& 

I 


k 


i 


$ 


$ 


I 

Jr* 


Rf 

§ 


w 

ft: 

8 


RHP 


Thus  S^(l)  and  Sq ( 1 )  are  regarded  as  (MSPE)^  and  (MSPE)^  in  the  notation 

2  2 

of  the  previous  section.  If  S^(l)  >  Sq(1),  we  then  say  that  the  level  0 
of  the  factor  has  less  contribution  to  MSPE  and  therefore  would  be 
preferred  to  the  level  1  in  view  of  stability  against  noise  factors. 

We  denote  Rank  R  Dj  R  ■  Vj  and  Rank  R  DQ  R  *  VQ.  We  now  present  the 
the  second  set  of  measures  of  dispersion  effects  of  levels  of  the  factor 
as 

S*(2)  -  (jr*  RD;R  , 

Sq(2)  -  (£'  R  DQ  R  £)/V0  ,  (5) 

at  the  levels  1  and  0,  respectively.  We  have 

MSE  -  (^)  Sj(2)  +  (^y)S*(2). 

2  2 

We  now  investigate  the  situation  where  S^l)  *  S^d),  u  “  0,1.  In 
other  words,  we  like  to  characterize  designs  for  which  y^  ■  y^.  We 
denote  the  row  of  the  matrix  X  corresponding  to  the  treatment  i  by 
x^(l*p)«  Note  that  for  each  i,  i-l,...,n,  the  row  is  repeated  r  times 
in  X.  Let  X*(n*p)  be  a  matrix  whose  ith  row  is  jc^.  Notice  that  rows  of 
X*  are  infact  distinct  rows  of  X.  We  have  X'X  -  r(X*'X*). 

Theorem  1.  For  i*l,...,n,  y^  «  y^  if  and  only  if  X*(X*'X*)  *  X*'  *  1^. 
The  proof  is  in  Appendix. 

Corollary.  For  n  ■  p,  we  have  Cl)  ■  S^(2),  u  ■  0,1. 

We  thus  observe  that  for  designs  with  n  ■  p,  two  sets  of  measures  are 
identical.  The  class  of  designs  with  n  ■  p  includes  the  known  Plackett 
and  Burman  designs  (see  Plackett  and  Burraan  1947).  We  however  strongly 


-6- 


6 


feel  chat  this  class  of  designs  is  very  weak  in  view  of  measuring 
dispersion  effects,  particularly  for  the  condition  that  there  is  no 
significant  lack  of  fit* 

3.2  Adjusted  Residuals  and  Third  Set  of  Measures. 

We  denote 


where  Is  the  vector  of  all  observations  corresponding  co  treatments 
with  6^  ■  u,  u  ■  0,1  for  the  chosen  factor.  Two  vectors  of  residuals 
1.  rQ  1.  at  the  levels  1  and  0  of  the  factor  are  generally  cor¬ 
related  under  che  model  (1-3).  We  now  present  a  vector  of  "adjusted 
residuals"  at  the  level  0  of  the  factor,  adjusted  w.r.t.  r^  so  that 
it  is  uncorrelated  with  We  denote 


/rll\  _  /Rlll  R1 12  \  _  /  R1 01  \ 

\r  12 /  V  R 1 12  R1 13/  VR102/ 

where  Rm  ( V ^ xVj }  with,  its  rank  Vj,  r^  (Vi*n)  with  its  rank  V^.  In 

» 

fact  we  have  Rm  *  tllrll*  We  now  write 


0a 


r0  _  R'l01 


R~l  r 
111  11* 


(6) 


It  can  be  seen  that  Rank  r  -  [(N-p)-V  ]  ■  V  (say)  and  furthermore, 

0a  1  0a 

Cov(r^  jr,  r^a  y)  ■  0.  We  call  r^  the  vector  of  "adjusted  residuals" 
at  the  level  0  of  the  factor,  adjusted  w.r.t.  the  residuals  at  the  level 
1  of  the  factor.  We  denote 


'  0a 


r0al 

vr0a2 
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where  r^j  (VQax*0  *ts  rank  Vga*  We  now  *lave  che  sum  s9uare8  of 

the  sets  of  linear  functions  r^  £  and  r^^  ^  [see  Scheffe  1959]  as 

“('nil  -£'rii  lriicul"lrn  z  • 

SSKal  m  i’rOal  tr0alr0al]’lro.l  i  •  (7) 

with  d.f.  and  V^,  respectively.  We  present  the  measures  of 
dispersion  and  adjusted  dispersion  effects  of  levels  of  the  factor 

Sj<3>  -  (ss(ru  i)/Vi]  , 

S^(3)  -  [SS(C041  jJ/vJ  ,  (8) 

at  the  levels  1  and  0  (adjusted  for  level  1),  respectively.  We  have 


•(mr)  st(3>  +te 


=?y)  soaC3)- 


Following  the  above  approach  we  find  the  vector  r  £  of  adjusted 

X  o 

residuals  at  the  level  l  of  the  factor  adjusted  w. r.t.  rQ  £  so  that  it  is 

uncorrelated  with  r_  y.  Let  r,  ,(v,  xn)  be  a  submatrix  of  r.  such  that 

0  lal  la  la 

Rank  r.  .  ■  Rank  r.  «  V.  ,  r_.fvnxN)  be  a  submatrix  of  rn  with  Rank  r,., 
lal  la  la  01v  0  '  0  01 

•  Rank  rQ  ■  Vq.  We  again  present  the  measures  of  dispersion  and  adjusted 
dispersion  effects  of  levels  of  the  factor 

SIa(3>  ■  tSS(rlal  *>/V  J  • 

SJ<3>  -  [sSCr0i  j-)/v0]  ,  (10) 

at  the  levels  1  (adjusted  for  the  level  0)  and  0,  respectively.  We  have 


—  ^  S2 
<-p) j  la 


(3)  +  (oRo)  s0(3)  *  (N-P>  -  Vla  +  v0  -  V1  +  V0a-(11) 


Theorem  2.  The  following  results  are  true. 

i.  v,  >(S  fiJU-l),  V  >  (  \  (1-6  ))(r-l)  , 

la  i-i  1  ua  i-i  1 
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ii*  VlaSla(3)  -  (  1 

To.so.(3)i^ti1<1-4t)5(r-1)so(1)  - 

iii.  If  Vla  -  (  z  <Si)(r-l)  then  S^(3)  -  S2(l)  , 
i-1 

n 

iv.  If  VQa  -  (  Z  (l-«1))(r-l)  then  S2fi(3)  =■  Sg(l). 

We  now  study  the  measures  in  two  extreme  situations:  (1)  r^  and 
Tq  2.  are  uncorrelated,  i.e.,  R^q  «  0,  (ii)  r^  £  and  r^  jr  are  completely 
correlated,  i.e.,  r^  **  D  r^  for  some  matrix  D. 

Theorem  3.  Consider  the  situation  R.„  “  0.  Then  S2(3)  =»  S2(2)  *  S2  (3), 
- -  10  u  u  ua 

u-0, 1 . 

Theorem  4.  If  r„,  *  D  r..  then  we  have  r„  **  0,  V„  ■  0  and  SS(r„  .y) 
-  01  11  0a  0a  0a  1^* 

-  0. 

Theorem  3  tells  that  in  case  Rjo  ■  0  there  is  no  need  for  the  adjustment 

of  residuals.  Theorem  4  tells  that  in  case  r^  £  is  linearly  dependent  on 

dependent  on  r^  £  then  the  level  1  of  the  factor  makes  all  contribution 

to  SSE  and  the  level  0  does  not  make  any  additional  contribution  to  SSE. 

In  case  V  *  V,  *  0,  we  have  V.  *  V,  =  (N-p),  r„,  =  D  r  and  D  is 
Ua  la  U  1  01  11 

nonsingular.  This  is  a  situation  where  the  levels  0  and  1  have  equal 

dispersion  effects  because  of  the  design  influence.  It  follows  from 

Theorem  2  that  for  r  >  1,  Vq  and  are  both  nonzero.  (We  assume 

naturally  that  there  is  at  least  one  6^  ■  1  and  at  least  one  (1-6^)*  1.) 

For  the  case  r  ■  1,  at  least  one  of  and  V,  could  be  zero  or  both  of 

0a  la 

them  could  be  nonzero.  We  now  consider  the  important  situation  when 
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both  V_  and  V,  are  nonzero.  If  S^(3)  >  Maximum  (s^(3),  S~  (3)),  then 
Ua  la  l  u  ua 

the  level  0  of  the  factor  has  an  advantage  edge  over  the  level  1.  On  the 
other  hand  if  S^(3)  <  Minimum(SQ(3),  S*  (3)),  the  level  1  of  the  factor 
is  superior  to  the  level  0  in  terms  of  smaller  dispersion  effects. 

Example  1 

We  consider  the  example  from  Box  and  Meyer  (1986),  page  20,  and 
Taguchi  and  Wu  (1985),  page  68.  Daniel's  normal  probability  plot 
indicates  that,  over  the  ranges  studied,  only  factors  B  and  C  affect 
tensile  location  by  amounts  not  readily  attributed  to  noise  (see  Box  and 
Meyer  1986).  We  now  fit  the  following  standard  linear  model  to  the  data 

E(y(xi,X2))  -  U  +  ai  B  +  a2  C, 

where  x^  ■  0,1,  *  (2x^  -  l),  y  is  the  general  mean,  B  and  C  are  the 

main  effects  of  the  factors.  We  can  write  the  above  model  in  the  form 

n  n 

(1-3).  Notice  that  N  =  16,  n  **  4,  p  »  3,  (  £  6  )  -  [  E  (1-6  ))  *  2  for 

i-1  1  i-1  1 

both  factors  B  and  C.  The  F  value  for  the  lack  of  fit  test  under  the 

assumption  of  normality  is  .1971  (<  1)  and  we  therefore  conclude  that 

there  is  no  significant  lack  of  fit.  It  can  be  seen  that  Rio  *  0  for 

both  factors.  We  observe  that  V,  *  ■  7,  V,  *  V.  *  6.  It  follows 

10  la  0a 

from  Theorem  2  that  S*  (3)  =  S*(l)  and  S*  (3)  -  S^(l).  In  Table  1  we 

la  i  ua  u 

display  numerical  values  of  various  measures  of  dispersion  effects  for 
factors  B  and  C.  We  find  that  the  level  1  has  a  lower  dispersion  effect 
than  the  level  0  for  both  factors.  We  also  observe  that  the  discrepancy 
between  dispersion  effects  at  the  levels  0  and  1  w.r.t.  all  measures  is 
more  for  the  factor  C  than  the  factor  B. 


i 
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Table  1 

Numerical  Values  of  Measures  Of  Dispersion  Effects  For  Factors  B  and  C 


» 


Factor 

S1<1>  ■  Sla(3) 

S0<‘>  •  s0.(3) 

S^(2) 

s0<2) 

S2  (  3 ) 

so(3) 

B 

.2863 

.3479 

.2498 

.3027 

B 

.3071 

C 

.0279 

.6063 

.0284 

.5241 

H 

.5286 

4 .  Properties 

We  now  state  some  properties  of  the  descriptive  measures  proposed  in 

Seection  3  under  the  model  (1-3).  We  first  observe  that  the  measures 
2  2 

S^l)  and  Sg(l)  do  not  depend  on  the  fitted  model  and  all  other  measures 

2  2 

depend  on  the  fitted  model.  The  measures  S^(l)  and  ( 1 )  are  always 

2  2 

uncorrelated  under  the  model  (1-3).  The  measures  S^(2)  and  Sq(2)  may 
however  be  correlated.  Two  sets  of  linear  functions  of  observations  Di  R 
jr  and  Dg  R  are  uncorrelated  if  and  only  if  Dj  R  Dg  ■  0.  Therefore  if  D 
R  Dj  -  Ogthen  S2 ( 2 )  ^nd  S2 ( 2 )  0re  uncorrelated.  We  have  the 
following  results: 

2 

Theorem  5.  Suppose  %  N(xj3,c  I).  A  necessary  and  sufficient  condition 
that 


(1) 

(2) 


£'  R  Di  R  Z  2 

- - - -  %  Central  X.z  with  d.f. 


X'  R  D0  R  X  2 

—  ■  — . -  Central  X  with  d.f. 


Trace  Rjj , 


Trace  Rgg, 


(3)  and  furthermore,  (1)  and  (2)  are  statistically  independent, 


Notice  that  Di  R  Do  *  0  if  and  only  if  Rio  -  0.  Moreover,  +  Vq  m 

(N-p)  if  Rio  ■  0  and  Vi  +  Vq  could  be  greater  than  (N-p)  if  Riq  *  0.  We 

2  2 

qu«_otion  the  use  of  estimators  S^(2)  and  Sq(2)  for  comparison  unless 
RlO  -  0.  We  of  course  realize  that  the  condition  Riq  *  0  is  too 
stringent  to  satisfy  even  for  one  out  of  m  factors. 


Theorem  6.  The  following  results  are  true. 


n  r 


n  r 


,E.  .£.5i(yi'fJrl)  '  ,£.  .‘.(‘-‘iX’ilT'll  *  “• 


i-1  j-1  J  i-1  j-1 

b.  If  for  the  factor  Rio  *  0,  then 


b.l.  Rn  and  Rqo  are  idempotent  matrices, 
b-2- 

b.3.  X'  R  -  0,  u  -  0,1, 
u  uu 


n  r 


n  r 


4  ,1  ,£.  siyi(hr^  * £,  ,£  (‘-‘ihMrh)  ■ 0 


i-1  j-1 


i-1  j-1 


2  2 

The  measures  S^O)  and  S(i_u)a^),  u  -  0,1,  are  always  uncorrelated 
The  reason  for  adjusting  residuals  is  to  obtain  uncorrelated  dispersion 
effects. 

Example  2.  We  consider  a  25  factorial  experiment,  i.e.,  m  -  5.  Let  the 
design  T(8x5)  be 


>i  Jt'  •»  >1' V J,1  u’  *»*  ««Vhf  ij  kVl^  l't  lVlV 


We  consider  a  situation  where  the  model  for  a  main  effect  plan  fits  the 
data  adequately  i.e.,  there  is  no  significant  lack  of  fit.  Therefore 
n  ■  8  and  p  ■  6.  We  note  that  it  is  sufficient  to  consider  the  distinct 
rows  of  X.  The  first  column  of  the  matrix  X*  has  all  entries  unity  and 
the  other  five  columns  of  X*  are  obtained  from  T  replacing  0  by  (-1).  We 


denote 


-2  -2 


-2  -2 


It  can  be  easily  seen  that 


18  -  X*(X*,X*)‘1X*'  -i  (~|£)  . 

It  now  follows  that  for  factors  1,  2,  4  and  5,  S?(2)  -  S^(2),  S*  (3)  « 

1  0  0a 

Sla^3)  "  °*  and  Sl(3)  "  So(3)  “  MSE*  For  the  factor  3,  R  -  0,  i.e, 
S^(2)  and  Sq(2)  are  uncorrelated  and  by  Theorem  3,  S^(3)  *  S^(2)  • 

S^a ( 3 ) ,  u  -  0,1. 

5.  Conclusions . 

In  industrial  experiments  for  quality  improvement,  dispersion 
effects  at  various  levels  of  factors  play  an  important  role.  They  are 
instrumental  in  the  choice  of  an  optimum  combination  of  levels  of  control 
factors.  This  article  presents  the  descriptive  methods  of  measuring  and 
comparing  dispersion  effects  at  the  preliminary  stage  of  investigation. 
The  outcome  of  such  comparisons  will  suggest  more  appropriate  complex 
models  for  further  investigation.  We  however  believe  that  the  imple¬ 
mentation  of  the  methods  discussed  in  this  article  will  result  in  highly 
informative  conclusions.  Although  in  this  paper  we  find  the  dispersion 


effects  of  levels  of  a  factor  using  one  factor  at  a  time,  the  same 
approach  can  be  used  to  find  the  dispersion  effects  of  level  combinations 
of  factors  using  many  factors  at  a  time.  Unless  the  number  of 
observations  is  sufficiently  large  in  every  cell,  the  reliability  of  the 
measures  using  many  factors  will  be  questionable. 


APPENDIX:  Proofs  Of  The  Theorems 


Theorem  1. 

A 

The  condition  y^  -  y^ ,  i-l,...,n,  holds  if 


(X'X)'1*! 


(i 

l  0  for  i  *  i ' ;  i,  i'e{l, . 


»nl 


The  above  condition  may  be  expressed  as  X*(X'X)  1  X*'  *  -j  1^,  or, 
equivalently,  X*(X*'X*)  *  X*'  ■  I  .  This  completes  the  proof. 


Theorem  2. 

A  A 

It  can  be  checked  that  Covfy,.,  y  -y  1  •  Covfy,,  y  -y  )  and 

ij  u  u  i  u  u 

therefore  Cov(yiJ-y1,  yu~yj  -  0.  Moreover,  Covfy^-y^  yuw~yu)  -  0, 
i  *  u.  It  now  follows  that  any  contrast  of  (y^-y^),  j  *  l,...,r  for 
a  fixed  i  with  6^  ■  1,  is  orthogonal  to  any  contrast  of  (yuw~yu)> 
w  ■  l,...,r  for  a  fixed  u  with  (l-6u)  *  1»  Furthermore,  any  contrast  of 
(y^-y^)  for  all  i  with  6^  -  1  la  orthogonal  to  any  contrast  of  (yuw~yu)» 
w  ■  l,...,r,  for  a  fixed  u  with  ( 1-d  )  ■  1.  The  results  (i-iv)  follow 
immediately  from  the  above  facts,  the  relationship  between  the  rank  and 
the  number  of  orthogonal  contrasts,  and  the  fact  that  the  sum  of  squares 
is  equal  to  the  sum  of  sums  of  squares  of  orthogonal  contrasts. 

Theorem  3. 

We  first  show  that  S^(3)  -  S^(2),  or,  in  other  words,  SS(r^  * 
jr'  R  Dj[  R  We  observe  that 


-1 
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-  r’ 


trnr'll) 


-1 


V 


-  z.'  Tj  where  R  z  -  rJ  £  *  ^  R  jr  , 

Rn  £  “  (Ru  £)’  (Ru  *) 

-  (Dj  »i)'  (dx  Ri)  -  i'  R  Dl  R  2- 

It  follows  from  (9)  and  (11)  that  VQ& sj^(3)  -  jr’  R  DQ  Rit  i.e.,  S*  (3)  - 
2 

Sq(3),  and  consequently  Vg  -  Vga .  The  rest  is  similar.  This  completes 
the  proof . 

Theorem  4. 

We  write  r^  ■  Dg  r ^  for  a  matrix  Dg  whose  independent  rows  are  rows 
of  D.  This  implies  that  -  DQ  and  thus  DQ  -  R”^.  Hence 

from  (6)  we  get  rQa  -  0.  The  rest  is  clear.  This  completes  the  proof. 


Theorem  5. 

It  is  known  (see  Rao  1973)  that  a  necessary  and  sufficient  condition 
of  (1),  (2)  and  (3)  to  be  true,  is  that  (r  r)  (r  Dq  r)  ■  R  R  Dg  R  * 
0.  Now 


We  see  that  R 
implies  Rjg  ■  0. 


R  Dj  R  Dg  R 

R  Dg  R  -  0 
Again  R1q 


R11 

Rio 

R01 

R11 

R10 

O 

O 

06  , 

R 

R 

R 

Rr 

R 

01 

10 

01 

01 

10 

00 

j 

implies  Rgj  R^g  Rg^  ■  0  and  this  in  turn 
■  0  implies  that  R  R  Dg  R  ■  0.  This 


completes  the  proof . 
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Theorem  6. 

A 

The  fact  X'R  ■  0  implies  X'R  jr  •  X'(^£.)  ■  0.  Considering  the 
columns  of  X  for  the  general  mean  and  the  factor  chosen,  the  result  (a) 
follows.  The  results  b.l  and  b.2  follow  directly  from  the  structure  of  R 
and  the  fact  that  R  is  an  idempotent  matrix.  The  result  b.3  follows  from 
X'R  -  0.  From  b.3,  we  get  jJ'X^  Ruu  -  0,  i.e.,  ^  Ruu  Zy  “  °* 

A  A 

result  b.2  implies  that  ■  0  and  hence  the  result  b.4  is  true. 

This  completes  the  proof. 


•wvvvs 
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